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Abstract. We study forward neutrino-induced meson production off the nucleon in the resonance region. Our calculation 
is based on a dynamical coupled-channels (DCC) model that reasonably describes n(y)N —¥ KN ,r\N ,KA,K~L data in the 
resonance region. We apply the PCAC hypothesis to the DCC model to relate the izN reaction amplitude to the forward 
neutrino reaction amplitude. In this way, we give a prediction for vN — > nN, nnN , rjN, KA, XX reaction cross sections. 
The predicted vN — > nnN , rjN, KA, XT cross sections are, for the first time, based on a model extensively tested by data. 
We compare our results with those from the Rein-Sehgal model that has been very often used in the existing Monte Carlo 
simulators for neutrino experiments. We find a significant difference between them. 
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INTRODUCTION 

A groundbreaking measurement of non-zero 613 opened a new avenue for neutrino oscillation experiments. Next 
generation experiments will be targeting the leptonic CP violation and the mass hierarchy of the neutrino. To achieve 
this goal, it is essential to understand the neutrino-nucleus interaction more precisely, 10% or even better, over a rather 
wide kinematical region that covers quasi-elastic, resonance, and deeply inelastic scattering (DIS) regions. 

In this contribution, we are concerned with the resonance region, from the A(1232) through 2nd and 3rd resonance 
regions, up to W < 2 GeV. Several models have been developed for the neutrino-induced single pion production off 
the nucleon in the resonance region, and have been used as basic ingredients to construct neutrino-nucleus interaction 
models. Some of them considered coherent sum of resonance contributions [1, 2, 3, 4]. Some others additionally 
took account of non-resonant mechanism of the tree level [5, 6, 7]. Two of the present authors further considered 
the rescattering also so that the nN unitarity is maintained [8, 9]. So far, most models deal with only the single pion 
production. However, the neutrino-nucleon interaction in the resonance region is a multi-channel reaction. Two-pion 
production has a comparable contribution to the single pion production. 77 and kaon productions can also happen. 
In order to deal with this kind of multi-channel reaction, an ideal approach is to develop a unitary coupled-channels 
model; this is what we will pursue. 

Recently we have been developing a unitary dynamical coupled-channels (DCC) model that can be extended to 
the neutrino reactions [10, 11]. Our DCC model is based on a comprehensive analysis of 71N, jN — > 71N, r]N,KA,KY. 
reactions in the resonance region, taking account of the coupled-channels unitarity including the %%N channel. In this 
contribution, we report our first step of extending the DCC model to the neutrino reaction [12]. For that, we invoke 
the Partially Conserved Axial Current (PCAC) hypothesis that allows us to relate cross sections of the pion-induced 
meson productions to those of the corresponding neutrino-induced meson productions at forward limit. 

In what follows, we give a brief description of our DCC model, and formulate the neutrino reaction cross section 
with the PCAC hypothesis. Then we show our result together with a comparison with the Rein-Sehgal model. 
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FIGURE 1. Comparison of the DCC model [10, 11] and data for single pion photoproduction. The upper [lower] figures are total 
cross sections (cr), unpolarized differential cross sections (do/dO.) and photon asymmetry (£) for yp —¥ 7t l) p [yp — > 7l + n\. The 
total energy is denoted by W, and the scattering angle of the pion by 6. 

DYNAMICAL COUPLED-CHANNELS MODEL 

In our DCC model [10, 11, 13], we consider 8 channels: yN, nN, T]N, 7lA,pN, aN,KA,KT,. The nicN channel is 
included in the nA,pN, ON channels using Feshbach's projection method, thus maintaining the three-body unitarity. 
Meson-exchange driving terms are derived from meson-baryon Lagrangian. The driving terms as well as bare N* 
excitation mechanisms are implemented in a coupled-channel Lippmann-Schwinger equation from which we obtain 
unitary reaction amplitudes. We analyzed n(y)N — > 7zN, r]N,KA,KY. reactions data simultaneously up to W = 2.1 GeV 
(W : total energy). The analysis includes fitting about 20,000 data points. As an example for showing the quality of 
the fit, we show in Fig. 1 the single pion photoproduction observables from the DCC model compared with data. As 
shown in the figure, our DCC model gives a reasonable description of meson production data in the resonance region. 
As a consequence, the DCC model contains all four-star resonances (and more) listed by the Particle Data Group [14]. 
Thus the DCC model provides a good basis with which we proceed to the neutrino reactions. 



PCAC AND NEUTRINO-INDUCED FORWARD MESON PRODUCTIONS 

Kinematic variables used in the following discussion are as follows. We consider the inclusive l(k) +N(p) —> 
l'(k') +X(p') reactions (X — nN,K%N,r\N... etc.), where (/,/') = (v e ,e~), (%,e + ) for the charged-current (CC) 
reactions. Although we do not show a result, the neutral-current reaction can be studied in a similar manner. We assume 
that leptons are massless. In the laboratory frame, the four-momentum are defined to be k = (E,k), p = (m^, 0,0,0), 
k' = (E' ,lc) and p' = k + p — k' . With the momentum transfer between I and I' , q = k! —k = (co,q), we define the 
positive quantity Q 2 by Q 2 = —q 2 = 4-EE' sin 2 j, where 9 is the scattering angle of I' with respect to I in the laboratory 
frame. 

For later use, we also define another frame where X is at rest. In this frame, q and p are denoted as q = (co c ,q c ) and 
p = (E^,—q c ), respectively, where = yj m 2 ^ + \q c \ 2 and m# is the nucleon mass. Also, we set q c = (0,0, \q c \) so 
that q c defines the z-direction of this frame. 



The cross sections for the inclusive neutrino and anti-neutrino reactions are expressed as 
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where the label a = CCv,CCv specifies the reactions; Q.' is the solid angle of /' in the laboratory frame; V ud is the 
CKM matrix element; the sign in front of Wj. a is taken to be + (— ) for v (v) induced reactions. The structure functions, 
Wi,a (i — 1,2,3), are Lorentz-invariant functions of two independent variables, (Q 2 ,W), where W is the total energy 
of X at its rest frame. In the forward limit, 0^0, Eq. (1) reduces to 
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The structure function W 2 . a is expressed in terms of matrix elements of weak currents between the initial nucleon 
N and the final state X, (X \j% \N) , as 
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where the summation symbol indicates all possible final states X, integration over momentum states of X, the average 
of initial nucleon spin state, and some kinematical factors including the phase-space factor. In the forward limit where 
Q 2 = 0, what survives in Eq. 3 is only the last term that contains the divergence of the current. The weak current 
consists of the vector (V) and axial (A' 1 ) currents. Because of the vector current conservation (X\q ■ V\N) = in the 
isospin limit, and \q c \ = (O c at Q 2 = 0. Thus we find 



W 2 , a (Q 2 ^0) = l i Y J \{X\ q -A\N)\ 2 . 

According to Refs. [15, 16, 17], we can define the pion field with the divergence of the axial currents as 

(X(p')\q.A a \N(p))=f Jt ml(X(p')\7i a \N(p)), 



(4) 



(5) 



where f n (%) is the pion decay constant (pion mass), and % a is the normalized interpolating pion field with the isospin 
state a. Furthermore, the matrix element {X{p')\jt a \N(p)) at Q 2 = can be expressed as 
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Here, ^7i a N^x(q 2 ) is the T-matrix element of the 7t a (q) +N(p) -^X(p') reaction in the kN center-of-mass frame, and 
the incoming pion is off-mass-shell q 2 = ^ m 2 . Using Eqs. (5) and (6), we have at Q 2 = 0, 
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where O n <>N^x is the total cross section of the on-shell % a + N —¥ X reactions, and we have used the relation 
■9~7i a N^rx{q 2 = 0) <~ .%[ a N^x(q 2 = m 2 ), which is a consequence from the PCAC hypothesis and q 2 = (O 2 . Converting 
the isospin basis to the charge state basis, we finally have 
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FIGURE 2. The structure function FiiQ 1 = 0) for the neutrino-induced meson productions from the DCC model. The solid (red), 
dashed (purple), dash-dotted (green), two-dotted dash (blue), and two-dash dotted (orange) curves are for the izN, ktzN, t\N, KA 
and K~L reactions, respectively. The sum of them is given by the thick solid (black) curve. The SL model is shown by the dotted 
(blue) curve. The figures are from Ref. [12], 



In the next section, we show the dimensionless structure function F2 defined by Fj_ = ©W2. The above result is 
essentially the same as those given by Adler [18] and also by Paschos et al. [19, 20, 21]. From Eqs. (2) and (8), 
one can evaluate neutrino-induced forward meson production cross sections at 9 = using the nN — > X total cross 
sections. 



RESULT 

We show the structure functions F% for the neutrino-induced meson productions off the nucleon in Fig. 2. The left panel 
of Fig. 2 shows CC neutrino-proton or antineutrino-neutron scattering where only 1 = 3/2 states (/ : total isospin of 
final state X) contribute, while the right panel is for CC neutrino-neutron or antineutrino-proton scattering for which 
both 1=1/2 and 3/2 states give contributions. While nN production is the dominant process up to W = 1.5 GeV, 
above that energy, the KizN production becomes comparable to nN, showing the importance of the nnN channel 
in the resonance region above A(1232). Also, we observe that the nN and nnN spectra above the A have rather 
bumpy structure, reflecting contributions from many nucleon resonances. This structure cannot be simulated by a 
naive extrapolation of the DIS model to the resonance region, as has been often done in previous analyses of neutrino 
oscillation experiments. Other meson productions, rjN, KA, and KL reactions have much smaller contribution, about 
[O(10 _1 )-O(10- 2 )] of 7TjVand nnN contributions. 

It is interesting to compare the above result for X = nN with the counterpart from the Sato-Lee (SL) model [8, 9], as 
done in Fig. 2. The SL model directly gives the Fi functions (without the PCAC hypothesis) because it consists of both 
the vector and axial currents, and reasonably reproduce available neutrino-induced pion production data in the A(1232) 
region [8]. From the comparison, we can see that contributions of nucleon resonances above A is clearly important 
above W = 13 GeV. Also, the good agreement W < 1.3 GeV indicates a reliability of calculating F2 from the nN — > X 
total cross sections with the PCAC hypothesis. We remark that this is the first prediction of the neutrino-induced nnN, 
r]N, KY production rates based on a model that has been extensively tested by data. 

It is also interesting to compare our result with F% from the Rein-Sehgal (RS) model [1,2] that has been used in 
many Monte Carlo simulators for analyzing neutrino experiments. Such a comparison is shown in Fig. 3. We can see 
that the RS model underestimates the A(1232) peak by ~ 20%. On the other hand, in higher energies, the RS model 
significantly overestimates Fi- Our result is based on the DCC model tested by lots of data in the resonance region 
while the RS model has not but based on a quark model. Considering that, the current Monte Carlo simulators using 
the RS model should be improved. In this work, the comparison with the RS model is done only at the forward limit. 
More comparison for non-forward kinematics, as well as full description of neutrino reaction needs development of a 
dynamical axial current model. Such a development is currently underway. 
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FIGURE 3. Comparison of F 2 (Q 2 = 0) between the DCC model and Rein-Sehgal (RS) model. 
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